Construction and analysis of a simplified many-body neutrino model 
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In dense neutrino systems, such as found in the early Universe, or near a supernova core, neutrino 
flavor evolution is affected by coherent neutrino-neutrino scattering. It has been recently suggested 
that many-particle quantum entanglement effects may play an essential role in these systems, poten- 
tially invalidating the traditional description in terms of a set of single-particle evolution equations. 
We model the neutrino system by a system of interacting spins, following an earlier work which 
showed that such a spin system can in some cases be solved exactly We extend this work by 
constructing an exact analytical solution to a more general spin system, including initial states with 
asymmetric spin distribution and, moreover, not necessarily aligned along the same axis. Our solu- 
tion exhibits a rich set of behaviors, including coherent oscillations and dephasing and a transition 
from the classical to quantum regimes. We argue that the classical evolution of the spin system 
captures the entire coherent behavior of the neutrino system, while the quantum effects in the spin 
system capture some, but not all, of the neutrino incoherent evolution. By comparing the spin and 
neutrino systems, we find no evidence for the violation of the accepted one-body description, though 
the argument involves some subtleties not appreciated before. The analysis in this paper may apply 
to other two-state systems beyond the neutrino field. 

PACS numbers: 



I. INTRODUCTION 

In many neutrino systems that are currently studied 
the rate of incoherent interactions is low enough to be 
completely negligible, yet coherent interactions (refrac- 
tion) play an important, even crucial role. A classical 
example is provided by the case of solar neutrinos: these 
neutrinos hardly scatter inside the Sun; nevertheless, 
their coherent interactions with the solar matter plays an 
essential role in their flavor evolution. This is of course 
the celebrated MSW theory of neutrino oscillations in 
matter HQ. 

The classical MSW theory describes neutrino propaga- 
tion in a background of normal matter (electrons, neu- 
trons, protons). There are systems, however, where the 
number density of neutrinos themselves exceeds those of 
electrons and baryons, such as the early universe, or the 
so-called hot bubble region in an exploding supernova. 
Such systems additionally require a theory describing 
neutrino self-refraction. 

Early investigations treated the neutrino background 
analogously to the MSW theory 0,0. Pantaleone @,T 
McKellar and Thomson [||, and Sigl and Raffelt 
showed, however, that the neutrino background is dis- 
tinguished by an important subtlety: the induced mass 
terms in general have nonzero off-diagonal components 



'Electronic address: 
T Electronic address: 
* Electronic address: 



fricdland@lanl.gov 

b.mckellar@physics.unimelb.edu.au 

ivona@physics.unimelb.edu.au 



in the flavor basis. Physically, this means that flavor can 
be coherently exchanged between the neutrinos. 

The authors of IE IS E| constructed the equation 
of motion of the neutrino system by using a one-body 
description for each neutrino. This treatment crucially 
depends on the assumption that the state of the system 
can be factorized into a product of one-particle states. 
If this is not the case and the wavefunctions of individ- 
ual neutrinos are entangled, a very different treatment 
may be required. A priori, it is not obvious that neu- 
trinos would not develop such entanglement, due to the 
off-diagonal induced mass terms. So the question arises 
whether this entanglement exists and has a substantial 
impact on the coherent flavor evolution of the neutrino 
ensemble. The answer could have a significant impact 
on the predictions for the supernova neutrino signal, the 
synthesis of heavy elements in the supernova, and possi- 
bly Big Bang Nucleosynthesis and cosmology. 

The question was recently examined by Friedland & 
Lunardini (F&L I 0, F&L II Q) and by Bell, Rawl- 
inson and Sawyer (BRS [Tl|). All three papers used a 
similar setup, with a Hamiltonian that was restricted to 
"forward scattering" only Q . While F&L I argued that 
the coherent part of the neutrino evolution should be de- 
scribed by the one-particle formalism, the BRS paper 
reached an opposite conclusion. 

To make their argument, BRS considered the evolu- 
tion of a system initially in the flavor eigenstates. With 
this choice, the one-particle formalism predicts no coher- 
ent flavor conversion and thus conversion on the coher- 
ent time scale in this system would be an indication of 
the breakdown of the one-particle description (presum- 
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ably through the formation of many-neutrino entangled 
states). The numerical calculation performed by BRS 
seemed to suggest the presence of such "fast" conversions, 
although the calculations involved a relatively small num- 
bers of neutrinos. 

F&L II pj subsequently solved the neutrino model in- 
troduced by BRS analytically for the special case of equal 
numbers of each neutrino species and equal strength in- 
teractions by mapping neutrino-neutrino interactions to 
spin- spin interactions. The solution in the limit of many 
particles exhibited the equilibration time that is precisely 
what would be expected from incoherent scattering. The 
analytical solution thus supported the one-particle de- 
scription of the system. 

One may wonder, however, if the initial state consid- 
ered in F&L II, namely, equal numbers of spins "up" 
and "down" was somehow special. Could entanglement 
appears with a more general setup? 

In this paper we will present a generalization of the 
many-body neutrino model introduced in F&L II, in the 
hope of understanding the quantum system better. We 
again consider a system of many neutrinos in which there 
are two flavor species, so that it maps to a system of 
interacting spins (thanks to the SU (2) symmetry of the 
problem). We generalize the model to initial states in 
which the species are not equally populated and where 
the initial states are not necessarily in flavor eigenstates. 
We show that the corresponding spin problem can still 
be solved exactly. The resulting solution exhibits a rich 
set of behaviors, as will be discussed in the following. 

Even generalized in this way, the model still involves 
several simplifications and it is important to spell these 
out. 

• The model keeps only forward scattering terms. 
Thus, while the model should correctly capture co- 
herent effects in a real neutrino system, conclusions 
about incoherent scattering effects in our model 
must be interpreted with care. 

• The momentum degrees of freedom are ignored. 
The effects of Fermi statistics are not included, that 
is the physical neutrino system is assumed to have 
a phase-space density much less than one. This 
is indeed satisfied everywhere outside the neutri- 
nosphere in a supernova. The case of the neutrinos 
in the early universe may be more subtle . 

• The interaction strength between any two neutri- 
nos is taken to be the same, ignoring the angular 



1 For a thermal neutrino distribution with zero chemical potential, 
the phase space density is given by (exp(_E/fcT) + 1) . This 
means that while the neutrinos on the tail (E ^> kT) are non- 
degenerate, those with E < kT are mildly degenerate. In the 
later regime, the present approximation may be inadequate. Of 
course, if the chemical potential is significant, the phase space 
density will be of order one in some regions. 



distribution of the neutrino momenta (see later). 
The model thus aims to describe the physical sit- 
uation in an isotropic neutrino gas and hence may 
or may not capture all effects that could arise as a 
result of very anisotropic momentum distributions, 
such as those suggested in 0, E3 • 

These limitations and assumptions should be kept in 
mind when relating the results obtained for the spin sys- 
tem to the behavior of a real neutrino gas. 

The problem of the flavor evolution in dense neu- 
trino systems continues to receive a significant amount 
of attention. In addition to the above mentioned pa- 
pers 111 fiol [III IT2I IT3I ] , the reader is referred to 
llflj for recent progress. 

II. SETUP AND GOALS 

A. The Hamiltonian and Eigenvalues 

We follow and Q by considering a system con- 
sisting of interacting massless neutrinos represented by 
plane waves in a box of volume V. Since our primary 
motivation is to investigate coherent effects in the neu- 
trino system, in particular, the possible breakdown of 
the one-body approximation due to flavor exchange, and 
not due to spatially dependent many-body correlations, 
we focus on the "forward" neutral current interactions 
between the neutrinos. In other words, we drop the mo- 
mentum degrees of freedom and include only scattering 
events that preserve neutrino momenta and those that 
exchange the momenta |g, [10| , 

\V x (k)Vy{p)) ► \V X {k)Vy{p)), (1) 

\vx{k)v y {p)) — > \v x (p)v y (k)) . (2) 

In the "usual" case of electrons, protons and neutrons 
in the background, the waves scattered forward interfere 
coherently. For the neutrino background, however, this is 
not necessarily so and our model, in addition to coherent 
effects, captures some of the incoherent effects as well 
PI ^j. The identification of coherent and incoherent 
effects will be discussed at length in what follows. 

For E v -C tti^o the neutral current interaction Hamil- 
tonian is 



Hint = [^Vort^Vaj ^y^V^Vpj . (3) 

Here the sum is over all k flavors. The Hamiltonian 
is invariant under a flavor SU (k) symmetry. Let us con- 
sider only two neutrino species, in which case the symme- 
try becomes SU(2) and the flavor space structure of the 
v — v interaction becomes equivalent to the interaction 
between pairs of spins. As explicitly shown in F&L II, 
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the interaction energy of two neutrinos, 1 and 2, is pro- 
portional to the square of the total angular momentum 
of the corresponding spin system, 



L 2 



(a 1 /2 + a 2 /2) 2 = (3/2 + ^-^/2). 



(4) 



The coefficient of proportionality dependents on the 
relative angle between the neutrino momenta 0, g = 
v^Gf ^ _ CO s©) Hj. i n a realistic neutrino system, the 
couplings are distributed according to the distribution of 
the relative angles between neutrino momenta. In order 
to make our system solvable, we will simplify the problem 
and take all the couplings to be the same 



9 = 



(•5) 



It is hoped that this simplification preserves the essential 



features of the evolution (see, however, 013). 

Henceforth we study this system of interacting spins to 
obtain information about the neutrino system. 

We will consider a system of N + M spins, such that 
initially N spins all have a certain orientation (for defi- 
niteness, without a loss of generality, "up") and the re- 
maining M spins all have a certain different orientation. 
At t = 0, thus, the N spins combine in a state of angular 
momentum Jn — N/2 and projection = N/2, and 
the M spins in a state of angular momenta Jm — M/2 
and projection hum- ln terms of the original neutrino 
system, we have a system of N electron neutrinos, \v e }i 
and M neutrinos in some other state \v x ). We give the 
answer for a general \v x ) and then explicitly study two 
cases: \v x ) is flavor eigenstate \v^) and \v x ) is a flavor 
superposition state (\v e ) + \uy))/\/2. 

The Hamiltonian for this system is 



J 



N+M-l N+M 



N+M-l N+M 



H = a E £fa + **) a= £ £ 2a r a k + -(N + M)(N + M-l) 



(6) 



j=i k=j+i 



j=i k=j+i 
I 



which is related to the square of the total angular mo- 
mentum of the system Q , 



' N+M 



j 2 = E ** 



0=1 

'N+M-l N+M 



E E 2a r a k + l (N + M)]. (7) 

3 = 1 k=j+l 



By comparing Eqs. © and we find 

Hint = 9[J 2 + -AM + N)(M + N-2)}, 



(8) 



with eigenvalues 

E(J, N, M) = g [,J(J + 1) + 3(M+A0(M +j v-2) j ; (Q) 

where 



state. As discussed in the time scale with which this 
probability evolves, t cq , tells us whether the evolution has 
coherent or incoherent nature. In particular, for M = N, 



(gN)- 1 

(gVNY 



(13) 
(14) 



for coherent and incoherent evolution correspondingly. 
In a large spin system the coherent time scale is much 
shorter then the incoherent one. One of our goals will be 
to see which timescales are present in our solution under 
different initial conditions. 

The second goal is to compare the coherent evolution 
we find to the predictions of the one-particle formalism. 
According to this formalism, the coherent evolution of a 
given neutrino is determined by the following one-particle 
Hamiltonian 0, 0: 



(10) 



J min = max.{\m N -m M \,\N/2-M/2\), (11) 
Jmax = N/2 + M/2. (12) 



Hi. 



part 



V2Gi 
V 



(15) 



B. Goals 

We are interested in finding the probability, as a func- 
tion of time, -Pi(t), of one of the particles remaining in 
the "spin up" state if it was initially in the "spin up" 



Here, ipi is the flavor state of the ith "background" neu- 
trino i.e. the background is all the neutrinos except for 
the one for which the equation is written. Explicitly, for 

T v l 1 

two neutrino species (y e and Vn) Wi) = f ■ The sum 
runs over all "background" neutrinos. 
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III. THE PROBABILITY OF SPIN 
PRESERVATION 

A. Result 

As we show in this Section, the evolution of our system 
can be solved exactly. In the interests of clarity, we begin 
by displaying the answer for Pi (t) : 

Pi® = \ + E C ( J ) 

J=Jmin 
./max-1 

+ £ r/(J) cos[0t(2J + 2)]. (16) 

J=Jmin 

I 



C(J) = 

x 
x 



ry(J) = 

x 

x 



The inner products in the last two equations are 
Clebsch-Gordan coefficients, the objects in the parenthe- 
ses are 3j-coefficients and those in curly brackets are 6j- 
coefficients. For definitions, see, e.g., |20| . 

The rest of this Section presents two complementary 
derivations of these results. The derivations are some- 
what technical and the reader primarily interested in the 
analysis of the rich physical properties of the solution 
may wish to skip to Sect. IIVI and return to this Section 
later, as needed. 



B. Construction of the probability: overview 

This solution can be found in either of two ways, both 
of which provide important, complementary physical in- 
sights into the spin system. These insights will prove 
very useful later, as we discuss the physical properties of 
the solution. Additionally, one or the other method may 
be useful for addressing still more general spin configu- 
rations. Correspondingly, we show both methods. 

• The first approach is to "split off" the first spin 
from the system, so that the remaining N + M — 1 
spins forms "a background" it interacts with. The 
solution is constructed by first coupling the angular 



The first two terms in Eq. i|16|l give the mean value 
of the probability: Pi — 1/2 means complete depolar- 
ization, or in the language of the neutrino system, an 
equal incoherent mixture of the two flavors; the second 
term thus given the degree of polarization of the mean 
("equilibrated") state. The last term contains the time 
evolution of the system. 



The limits of the summation J m i n and J max are given 
in Eq. {TTJ and (T2J. The coefficients C(J) and T](J) 
given by 



(17) 



(18) 

I 

momenta of the remaining spins and then coupling 
the first spin to the result. This method generalizes 
the idea employed in 1] , but without relying on the 
symmetries specific to the N — M case. 

• The second approach is to treat the first spin as 
a part of the system. We solve for the evolution 
of all spins that start out in the "up" state. The 
solution can be found by observing that even with 
N M the system possesses a very high degree 
of symmetry: all spins that start out in the same 
state evolve in the same way. This means, as we 
will show, that the problem can be reduced to that 
of just two coupled angular momenta. 



C. Method I: splitting off the first spin 

The outline of this approach is as follows. The time 
evolution of the system is easily written down in the ba- 
sis of total angular momentum, since, as explained in 
Sect. 1^ this is the Hamiltonian eigenbasis. Correspond- 
ingly, we begin in Sect. lIII C ll bv constructing the density 
matrix for the system in this basis. In Sect. IIII C 21 this 
density matrix is rotated to a basis in which the first 
spin has a well-defined value. Finally, in Sect. IIII C 31 the 



(_ 1 )2J-( mfl+ ™„)+(jv+M)/2 + ij 1 + 2 J) — y/(N/2)(N/2 + 1)(N + 1) 

\(N/2, m N , M/2, m M \J, m N + m M )\ 2 

( J 1 J \ | N/2 J M/2 1 

\-{m N + m M ) m N +m M )\ J N/2 1 /' 

^2J-(m N +rn M ) + (N+M)/2^^ + 2J ^g + 2J)^y/(N/2)(N/2 + 1)(N + 1) 

(N/2,m N ,M/2,m M \J,m N + m M ) (N/2, m N , M/2, m M \ J + ^,m N + m M ) 
( J+l 1 J \f^/2J+l M/2 1 
(mN + mM ) mN+TnMj\ J N/2 1 J' 
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probability P\{t) is found. 

1. Constructing the many-body density matrix in the total 
angular momentum basis 

Our system begins (at time t = 0) in the state 

S(0) = | j N , m N ) <g) | j M , m M ). (19) 

Recall jjv = N/2 is the total angular momentum of all the 
spin up particles, each with angular momentum i, and 
projection along the z direction, mjq. Also, ]m = M/2 is 
the angular momentum of all the background particles, 
each with angular momentu m ^ and projection, vhm- 
Rotating the initial state, Ea. 1)19(1 . to the total angular 



momentum (J) basis and evolving it to time t, we have, 
S(t) = ^2 (jN, hi ,J,m\j N , j M , m N , m M ) 

J=Jmin 

xe- itE ( J ' N M\j N ,j M ,J,m). (20) 

Here ( j Nl j M ,J,m \ jN,jM,m N ,m M ) is the Clebsch- 
Gordan coefficient where jn and jm are coupled to the 
total angular momentum of the system J with projection 
in the z direction of m = tojv + Tom- 
The density matrix is defined as 

p(t) = \S(t)){S(t)\. (21) 

Hence we have, 



P(t) = Ej™ j min E j"=J min e ltAE{J ' J ) ( ]N,jM, J, m | j N ,j M ,m N , tom ) ( 3n,3m, J', m \ j N ,jM, m N , m M ) 

\ jN,3M,J,m)( j N ,j M ,J',m |. (22) 




FIG. 1: Graphical representation of recoupling. The diagram 
on the left represents the original coupling: the net angular 
momentum (ji) of the iV spins pointing up and the net angu- 
lar momentum (j'2) of the remaining M spins are coupled to 
the total angular momentum of the system (J). The diagram 
on the right represents the new coupling (see main text for 
description), k is the momentum of — 1 "up" spins. 



Here AE(J,J') = g[J(J+l)-J'(J' + l)] is the difference 
between the eigenvalues E(J,N, M) and E(J',N,M). 



2. Construction of the density matrix in the basis 
1 1, A)® 

To find the probability of spin preservation it is conve- 
nient to again change the basis, this time to | h , A) ® \j, /x) . 
|i,A) is the state of one of the spin-1/2 particles with 
projection A and \ j, fi) is the state of M + N — 1 remain- 



ing particles with angular momenta j and projection /z. 
The density matrix in the new basis is constructed in 
this subsection and then the probability is found in the 
next subsection. We start with the state of Eq. (|2()|l , and 
transform it to the new basis. We will omit the limits on 
the summation signs for the next few equations but will 
comment on these later. 

To change to the preferred basis (|^,A) <g> we 
modify how the angular momenta couple to the total an- 
gular momentum, J. This can be done in any way that is 
convenient. We for the moment "remove" the first spin 
up and combine the N — 1 spins up (with net angular 
momentum k = jjy — 1/2) and the remaining M spins 
(with net angular momentum 3m) into an object with 
angular momentum 3. We then couple the removed spin 
and 3 to the total J. Therefore 3 is the result of coupling 
3m and k. See Fig. ^ for a pictorial interpretation. 

The state | 3N,3M,J,ni ) is represented in the new 
basis as 

I 3n,3m, J,m ) = ^( -, (k j M ) j; J \ (- k)j N ,j M ; J ) 
j 

x \\,(kj M )j;J,m). (23) 

The notation (a b)c indicates that a and b 

couple to c. Here the recoupling coefficient 

( \,{k 3m)3\J I (5 k)j N ,j M ;J ) is proportional to 
the 6-j coefficient, 
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< i, (k Jm )j; J | (i fc) jN , JM ; J ) = (_i)i+*+J+iM [(2jjy + 1)(2j + ^ JA,}. (24 ) 



For a definition and explanation of 6-j coefficients see 
|2lj or HJ, Eq. (108.6). The right side of this coefficient 
represents the original coupling (the left side of Fig.^l. 
The left side represents the changed coupling (the right 
side of Fig. QJ. The 6-j coefficient is a consequence of 



Finally, substituting Eq. into Eq. l(2U)l yields 



Recall that the first N spins are initially in a state with 
total angular momentum jjf = N/2. The initial state of 
the test spin is consequently the same as that of the other 
N - 1 particles, "up" (A = +1/2). 



I 

changing the coupling. 

Rotating from the basis | ju)j>J) m ) to the 

basis | \, A ) <E> | j, p ) and substituting k = jn — 1/2, we 
have, 



(25) 



r 



(26) 



It is now simple to construct the density matrix in this 
new basis using Eq. {gj and p = | S'(i) )( |. We 
present the density matrix in component form, 



I 

| j N ,j M , J,m) = J2 (-l) J+jM+j » [(2j N + l)(2j + 1)]' ( \, j, A,/i | 1 j; J,m ) 



I 

S(t) = £ (-l) J+J ' M+jV [(2iJV + l)(2j + l)] l e-«W' M )(i JV ,^ 

I 



Pi A,i„ ; A A' , ^ (*) = ^e~ ltAB(J ' J,) (-l) J " J '(2iJV + l)(2j + 1)( i2v,iM, J, m | jN,3M,m N , m M ) 

x (jN,3M,J',m I jN,3M,m N ,m M )( A,/i | -, j; J, m ) (27) 



1 1 



i i „• i r i i 



' V „' I _ A'. T> ™ \ J 2 — o, 3N I J 2 JN- 2 JN 



v 2' J ' ' p 1 2' 7 ' 7 1 JM J J J 1 JM J J 



3. Probability of Spin Preservation 

The probability that the first spin remains in the up 
state can be found from the density matrix. In general 
the probability of an eigenvalue a, represented by the op- 



erator A is Tr(Ap). In this case the operator is diagonal 
and hence probabilities are just the diagonal components 
of the density matrix. These components of the density 
matrix which give probabilities are those where j = j', 
p = pi ', and A = A', i.e., pi x j jU . i Aj j |U . Furthermore we 
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are looking for the probability of the first spin remaining 
in the spin up state so that A = +A. Using this infor- 

I 



mation together with Eq. (|27|l we find the probability to 
be 



Pi(t) = ^> A. , - 



E 



-itAE(J,j')(_ -\\J-J 



J> J,J':3,M 



(-1)' (2jjv + l)(2j + i)(jN,3M,J,m | 3N,jM,m N ,rnM ) 



x (jN,3M,J',m | 3n,3m,™>n,™>m )( -, j, ^ /x | J,m)( -,// | -,j;J',m) 

JM ■/ j J I J' j 



(28) 



The third Clebsch-Gordan coefficient in Eq. 
( | ^,j;J,m ), represents the coupling of | 

and j to J. This gives J— In 
this same Clebsch-Gordan notice that we must have 
A + /x = rn, so that the summation over /i is unneces- 
sary. Analyzing the fourth Clebsch-Gordan coefficient, 
( \.,3, 5,// | Tj,j;J',m ), further simplifies the equa- 
tion. This coefficient shows that we couple \ and J to J' . 
Hence we have, |j — Aj < J' + Now if j = J+| then 
J < J' <; J + f , and if j = J - \ then | J - 1| < J' < J. 
Therefore J' = \ J - 1\, J, J + 1. 



Summing over J' we find 



and 

7 ? (J) = -2(2j JV + l)(2J + 2) 

x ( 3n , Jm ,J,m\j N ,jM,m N ,m M ) 

x (3n,3m,J + 1,to I jN,3M,m N ,m M ) 



1 x 1 1 

2' J+ 2'2 / ' 

1 x 1 1 

2' J+ 2'2 / ' 

\ 3N~\ 
3m J 



1 

n 

2 

1 

n 

2 

3N 

J+i 



J,m' 



J + 1 , TO 



2 

3m 



JN 
J- 



JN 
J+i 



(31) 



Although it is not immediately obvious, the answer just 
obtained coincides with Eqs. lfffi|l . l(T7|. and l(T%)l. Notice 
that (7(7) is not the same as C(J) in Eq. (|f 7|l since, 
unlike Eq. (| 1 6fl , Eq. I|29|) does not contain the 1 /2 term. 



p 1 (t)= 

+ n(J) cos[tAE(J + 1, J)], 



where 



C(J)= (2jjv + l)(2j + 1) 

x I < JnJm, J, m I j N ,3M, m N , m M ) | 2 



( 2'-?' 2' m ~ 2 ' 2 ,J ' ' m ' 



5 JJV - I jjv 



(30) 



D. Method II: symmetry of the entire system 

The second method is based on an observation that 
since the coupling strength is the same for all spins the 
system possesses a very high degree of symmetry. In par- 
ticular, all spins that start out in the same state evolve 
in the same way (up to relabeling). This simple obser- 
vation proves to be very powerful for our analysis: it 
means that all spins that start out in the "up" state com- 
bine in a single composite object with the angular mo- 
mentum L\ — N/2 at any moment in time, not just at 
t = 0. The same can be said about the remaining M 
spins (L 2 — M/2). Thus the problem reduces to a text- 
book case of just two interacting angular momenta. 

More precisely, at any moment in time, the system 
is described as a superposition of states i with definite 
values of total angular momentum, Ji. Each such state 
is obtained as a result of symmctrizations and antisym- 
metrizations of individual spin wavefunctions. The high- 
est J state is a completely symmetrical combination, 
the next one is obtained by performing one antisym- 
metrization, etc (the general rules are given by Young's 
tableaux). The important point is that by symmetry 
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all spins that started out in the same state can only be 
symmetrized. That means that all spins that started 
out in the same state can be assembled in a single object 
with the total (unchanging in time!) angular momentum. 
This argument holds for any t. 

It is conceptually straightforward to compute the ex- 
pectation value of the operator jj^ (the z component of 
the angular momentum of the first composite object) at 
any moment in time. The expectation value of the corre- 
sponding operator for each spin, can then be found 
by simply dividing by N, by symmetry. 

Let us denote the state of the whole system by \S(t)). 
As before, the system evolves according to 



where J m in and J ma x are given in Eqs. (fTTfl and 
Ej = E(J, N, M)) given in Eq. © (the irrelevant con- 
stant piece can be dropped). 

We need to find the probability of the first spin being 
in the "up" state, P x (t) = (5(t)| +l/2) 1 (+l/2| 1 |5(i)). It 

follows from the definition of s [ P = 1/2(| + l/2)(+l/2| - 
| — 1/2) (— 1/2 1) and the completeness relation (| + 
l/2)(+l/2| + | - l/2)(-l/2| = 1) that this probability 
is related to the expectation value of the angular mo- 
mentum of that spin, 



P x {t) = (S{t)\s^\S{t)) +1/2. 



(33) 



\S{t))= J2 

•/=./mi 



\J){J\S(0))t 



-itEj 



(32) Let us now compute (S^s^l-SC*)): 



1 1 ;„« 



(34) 



./=,/„- 



J'=Jn 



The products (J\S(0)), (5(0) | J') are just the Clebsch- 
Gordan coefficients (N/2, mjy, M/2, mM, \ J, rnj^ + tom), 
{N/2, ttin, M/2, mu-, \J', wat + TOm) (written so because 
they are real). The problem of finding the expecta- 
tion value (J'\j^\J) is solved, e.g., in [2(]]. Using Eqs. 
(109.2), (109.3), (29.13) there, and the fact that the op- 
erator ji 1 ' is a q = component of a k = 1 spherical 
tensor (see (107.1) of |20j for the exact definition), we 
get 



J' 1 J 

~(m N + m M ) m N + m M 



The 6j-symbol enforces the following selection rule: 

J' = J,J±1. (37) 
The double sum in Eq. (|34|l then reduces to 

(S(t)\s^\S(t)} = 
1 

— \(N/2,m N ,M/2,m M ,\J,m N + m M )\ 2 



J=Jmi„ 

J 1 J 

~(m N +m M ) m N + m M 



x M + + { »? ^ M { 2 } * " + "»{T//2T } 



x y/(N/2)(N/2 + l)(N + l) 

x (-iy, 



(35) 



where 



x v / {N/2)(N/2 + l){N + 1) 

x ^_^2J-(m N +m N )+N/2+M/2+l 

— Y (N/2,m N ,M/2,m M ,\J, m N + m M ) 



N 



q = max[J, J'} — (m N + m N ) + N/2 + M/2 
+ min[J, J'] + 1. (36) 

The first line in Eq. I|35|) expresses the dependence of 
the matrix element on the ^-component of the angular 
momentum of the whole system (Eq. (109.2) of |20|). 
The second line contains the dependence of the reduced 
matrix element on the total angular momenta J and J' 
(Eq. (109.3) of 20]). The third line is the reduced matrix 
element of the angular momentum of the first composite 
object, (j«|j«|j«) - y /L 1 (L 1 + l)(l + 2L 1 ). Finally, 
the last line contains the sign factor collected from all the 
ingredients. 



J=J„ 



x (N/2, m N , M/2, m M , \ J + 1, m N + m M ) 
J + 1 1 J 
-(m N +m M ) m N + m M/ 

x #Tzlmi) { T ^ T } 

x y/(N/2)(N/2 + l)(N + 1) 

x ^_^2J-(m N +m N )+N/2+M/2 

x 2cos[2fift(J+l)], 



(38) 



which is the answer given earlier in Eqs. 110|) . (|17|l . and 
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IV. ANALYSIS: A FLAVOR SUPERPOSITION 
BACKGROUND 

As already mentioned, the solution in Eqs. ifTBj l. (jT7|l . 

and l|18fl contains very rich physics found in the actual 
neutrino system. First, let us show that it contains 
both coherent oscillations and incoherent equilibration 
(decay). For that, let us consider the case 



N 

m N = —, 



m M = 0. 



(39) 



In the neutrino system, this corresponds to the first 
N neutrinos initially being in the \v e ) flavor eigenstate, 
while the remaining M neutrinos starting out in the max- 
imally mixed state. In this setup, the one-particle formal- 
ism predicts coherent evolution. 



0.5 




200 



400 

x 



600 



800 



FIG. 2: Time evolution of Pi(t) in the case N = 100, M = 
120, mjv = N/2, rriM = 0. The time variable is defined as 
r = gt(M + N). 

Indeed, our solution confirms this. The explicit form of 
the coefficients C(J) and r]( J) is given in Eqs. (|A1() and 
(|A2|) . Using these values, we plot the behavior of the 
expectation value Pi(t) for N = 100, M = 120 in Fig. El 
We see that the evolution exhibits both oscillations and 
decay. The oscillations reflect coherent behavior, while 
the decay occurs on a longer time scale, as would be 
expected for incoherent evolution. 

This assertion can be made quantitative. The key test 
of coherence/incoherence of the evolution is the scaling 
of the evolution times with the number of particles. In 
the limit of large N and M, it is possible to show that 
a single "wavetrain" in Fig. [21 (oscillations plus decay) is 
described analytically by the following formula: 



1 / N 2 
1 / N 2 

+ 2 1 1 -FT^ |exp 



2,2 



gH 



4- 2-^ 



X cos 



(40) 



A detailed derivation of Eq. i|40|) is given in Appendix 
IB1 Here, we notice that the oscillation period scales as 
oc (N 2 + M 2 )~ x / 2 , as expected for coherent evolution, 
while the decay time goes like N~ 1 / 2 {2 + 2N 2 /M 2 ) 1 / 2 , 
indicating its incoherent nature (c/Eqs. I|13fl . CH)- 

We further observe that the frequency of the coherent 
oscillations agrees with the predictions of the one-particle 
formalism for a neutrino in the background of M maxi- 
mally neutrinos (\u e ) + \v^))/^/2 and N — 1 ~ N neutri- 
nos in the \v e ) state. Indeed, the one-particle oscillation 
Hamiltonian l|15|l in this case is 



Hi- 



■part — 



V2Gi 
V 



N - 



K 

2 



(41) 



so that the oscillation frequency, given by the difference 
of the eigenvalues, is precisely g^/N 2 + M 2 . 

A very valuable physical insight can be gained from the 
idea that underlies Method II of solving for the evolution 
of the system, namely, that all spins in the same initial 
state always combine to form an object with a certain 
definite value of the angular momentum. As explained 
in Sect. IjHI D|l . this means that the system can be re- 
duced to just two coupled angular momenta, ]n = N/2 
and 2m — M/2. When the numbers of spins N and 
M are sufficiently large, the two composite angular mo- 
menta behave as nearly classical objects. They precess 
about the direction of the total angular momentum of 
the system. These are the fast oscillations seen in Fig. |5] 

Quantum-mechanically, a system that has a definite 
value of the total angular momentum does not simultane- 
ously have definite projections of the individual angular 
momentum vectors that comprise it. Correspondingly, 
after a while the components of the angular momenta 
that are transverse to the total angular momentum un- 
dergo "quantum wash-out" . The system equilibrates to 
a state in which only the components of the constituent 
angular momenta along the direction of the total angular 
momentum remain. 

Let us check this quantitatively. The total angular 
momentum of the system (in the classical limit) makes 
an angle cos a = N / \/N 2 + M 2 with the positive z direc- 
tion. Hence, the projection of the net angular momentum 
of the first iV spins on the direction of the total angular 
momentum of the whole system is N/2 x cos a. After 
a sufficient amount of time this is the only component 
that remains, the transverse components are washed-out. 
Projecting it back on the z-axis, we get N/2 x cos 2 a or, 
using Eqs. pi jl . 



1 INN 2 

p, = _ -| 

2 N 2 N 2 + M 2 



N 2 



N 2 + M 2 



(42) 



precisely in agreement with Eq. Q40[l. 

Additional insight about classical and quantum fea- 
tures of the evolution can be gained by restoring the fac- 
tors of the Plank's constant in Eq. I)40[l . It is simple 
to see that the product gt comes with a one factor of 
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h. The logic is as follows: (i) from Eq. it is obvious 
that the energy is proportional to H 2 (from the angular 
momentum squared factor); (ii) next, in computing the 
evolution phase (exp(iEt/H)), we divide by one power of 
S; (iii) this leads to one power of h in the argument of 
the cosine in Eq. (14), i.e., gt comes with one factor of 
h. Eq. then reads 



N 2 
N 2 + M 2 

N 2 
N 2 + M 2 



exp 



-g 2 ht' 



Mh 

o (M i N 
\ N M ) _ 



X COS 



2 + M 2 t 



(43) 



We see that the argument of the cosine 
gh\/N 2 + M 2 t — gJtott involves only the classical 
value of the total angular momentum Jtot, while the 
decay exponent involves two powers of h, of which 
only one is absorbed into the definition of classical 
angular momenta. The decay exponent is a quantum 
effect, in a sense that its physical origin lies in the 
quantum uncertainty principle. This allows us to make 
an important identification: coherent evolution in the 
neutrino system maps into the classical behavior of 
the angular momenta, while incoherent effects in the 
neutrino system correspond to quantum effects in the 
spin system. This important point will be developed 
further in Sect. IV Al 

Finally, we note that the solution is periodic and the 
wavetrain, having seemingly completely decayed away, 
reemerges after some time. This phenomenon is due to 
the fact that the spin system possesses a fundamental 
frequency, of which all other frequencies (the arguments 
of the cosines in Eq. Ultjfl ) are multiples. This effect will 
be discussed in detail in Sect. IV Cl 



up (N) particles and a subset of these is shown in Fig. 01 
The time on the figures is scaled such that r = tg(N+M) 
so that we may compare our results to those in F&L II 
and BRS. 

The main features of the solution we find are the fol- 
lowing: 

• As seen in the left panel of Fig. for N ~ M (to 
be quantified later), the system equilibrates to some 
value of P\ and subsequently for a long time shows 
no evolution. This behavior is similar to what was 
observed in the case of Sect. (|IVI) . Just like that 
system, the wavetrain reemerges after some time. 

• As seen in the right panel of Fig. [3 the equilibration 
happens on time scales characteristic of incoherent 
evolution, t eq ~ £r x ((M + N)/2)~ 1 / 2 . Unlike the 
case of Sect. (jlVjl though, in the present case for 
N = M we do not see coherent oscillations with 
the decaying wavetrain. In fact, it was explicitly 
shown in F&L II that for TV = M the probability 
Pi depends on time only through the combination 
y/Wgt. Thus, the evolution has manifestly incoher- 
ent nature. (c/Eq. (fTOjl ). 

• As | N — M\ is increased, the value of P\ in equilib- 
rium, Pi, increases and for large \N— M\ the system 
stops evolving, entering a "freeze-out" state. This 
trend can be clearly seen in Fig. [3] As further illus- 
trated in Fig. 21 where Pi is plotted as a function 
of N — M for fixed M + N, the minimum value 
of Pi = 1/2 occurs when N = M. That Pi rises 
for both signs of (N — M) appears counterintuitive 
at first. Indeed, it means if we have a single spin 
up coupled to a very large system of spins down, 
the first spin does not equilibrate to a state mostly 
oriented down, but remains aligned up. 



V. ANALYSIS: A FLAVOR EIGENSTATE 
BACKGROUND 

We next analyze the probability for the case that the 
state jm is made of spin down particles (or neutrinos 
in the muon flavor eigenstate). Hence we identify the 
quantum numbers as 



N 

m N = y, 



m M = ~ — 



M 
2 



(44) 



One reason this case is interesting for us is that for this 
system the one-particle formalism predicts no coherent 
flavor conversion. Indeed, the off-diagonal terms in the 
Hamiltonian H15|l vanish. 

The explicit form of the coefficients C{ J) and rj(J) is 
given in Appendix |XJ Eqs. (|X3l) . JX5j|, and (|X4|) . As 
explained there, for the case N — M our probability 
agrees exactly with that found in F&L II. 

We have plotted the probability according to these 
equations for various numbers of spin down (M) and spin 



• As | TV — M\ is increased, in addition to decay, the 
system start exhibiting oscillations, the frequency 
of which grows with \N — M\. The amplitude 
of these oscillations gets progressively smaller as 
\M — N\ increases, since the system approaches 
the freeze-out state. This behavior can be seen in 
Fig. which reproduces at higher resolution the 
case M = 1280, N = 2600 of Fig. El 

In what follows, we will discuss these features further. 
This will allow us to gain a deeper understanding of the 
spin system and its relationship to the neutrino system. 
We will first show how the transition from the regime in 
which the spins equilibrate to the freeze-out regime illus- 
trates the interplay of "classical" and "quantum" effects. 
We will then discuss the time scale of the evolution and 
whether the system follows the predictions of the one- 
particle Hamiltonian in the classical regime. Finally, in 
Sect. IVTJl we will comment on the periodicity of our so- 
lution. 
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FIG. 3: Plots of N=2560 and various numbers of spin down particles, M. The time variable r is defined as the scaled time 
r = gt(M + TV). Note that for all the graphs N + M — even. The right figure shows that the equilibration time for the bottom 
two curves is characteristic of incoherent evolution, t eq oc [2/(M + N)} 1 ^ 2 . 
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FIG. 4: The average probability (Pi = 1/2 + ][]./ C( J)) when 
we change the number of spin ups and spin downs while keep- 
ing M + N constant. 



FIG. 5: A close up of the plot for N=2560 and M=1280. 
Note the vertical scale. The probability is very close to one 
and fluctuates minimally about the average value. 



A. From equilibration to freeze-out: a 
quantum-to-classical transition 

In the course of the analysis in Sect. IIVI we have en- 
countered a situation in which fast coherent processes in 
the neutrino system corresponded to classical effects in 
the spin system, while slower incoherent processes in the 
neutrino system corresponded to quantum effects in the 
spin system. This physical picture is further illustrated 
by the system — N/2, tum — —M/2 we are now 
considering. 

Once again, we recall that this system is reduced to 
two angular momenta as described in Sect. IIIIDl When 
N and M are large, these two angular momenta become 
approximately classical objects, ji and j<x- These objects 
are coupled with an interaction const x ji ■ j 2 - By en- 
ergy conservation, this quantity stays fixed, and since the 



lengths of both vectors are fixed the angle between the 
vectors j\ and j% stays fixed. Since the two vectors start 
out pointing in opposite directions, they remain point- 
ing in opposite directions. The only possibility left is for 
the system to tilt as a whole, but that would violate mo- 
mentum conservation, unless the total momentum of the 
system is zero (N — M). Thus, the classical system is 
frozen, unless ji — j%, which is a special point. 

Our calculation shows that quantum mechanical ef- 
fects actually resolve this discontinuity. The transition 
from the freeze-out to equilibration happens in some fi- 
nite range of angular momenta. 

What physics sets this range? Let us consider the two 
angular momenta being added as semiclassical objects. 
Then, in addition to the angular momenta in the z direc- 
tion, they each possess "quantum" angular momenta in 
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the x — y plane. These momenta are of the order 
J xy ~ V J 2 - Jf = y/J(J+l)-P = Vj. 
That this is a quantum effect is seen by restoring the 



units: J xy ~ v J/i. This "quantum" angular momentum 
is what is responsible for the equilibration. 

When the two angular momenta are added, their x — y 
quantum momenta are also combined and the net ob- 
ject has J x - y of the same order as the ingredients. If 
the net classical momentum along the z axis is greater 
then J x -y, the equilibration will be suppressed. In other 
words, we arrive at the physical condition that deter- 
mines the boundary between equilibration and "freeze- 
out" 



\M - N\/2 ~ y/(M + N)/2. 



(45) 



For large number of spins, we can simply replace in this 
condition sj (M + N)/2 -> VN. 

This boundary between freeze-out and equilibration 
can be quantified by the analytical expression for the 



equilibrium Pi that can be derived for \M 
(M + N) (see Appendix 0: 



Pi~l 



1 



F ( — 



Here, A = (N - M)/2, E = (M + N)/2 and 

i a 2 r°° p- y 

F(a) = - ~ - / dy 



y + a 1 



N\ < 



(46) 



(47) 



Up to a small correction in the prefactor in Eq. (|46J) . 
Pi depends on A and £ in combination A 2 /£. Moreover, 
P(a) quickly drops to zero as a is increased beyond one. 
Hence, the width \M — N\ of the "quantum equilibration 
region" indeed scales with TV + M as y/N + M. 

This scaling can be seen in Fig. 0] The curves shown 
there are obtained by using the exact solution, Eqs. (|16[) . 
QA3|) . \AA\. I|A5|) . They are well described by Eq. |@BJ|. 



B. Equilibration Time 

We now discuss the time scales that control the evolu- 
tion of our system. We are particularly interested to see 
how the evolution scales with the number of spins in the 
limit when the spin system is large. 

As already mentioned, the case of N = M considered 
in F&L II is straightforward: for large N the evolution 
is uniquely dependent on yN t and the decay curve is 

Pt(t) = 1 - ^VNgtexp{~Ng 2 t 2 )erfL{VNgt), (48) 

where erfi(z) is the imaginary error function 

2 



erfi(z) = —i erf(iz) 



x / exp(r)df . (49) 



Our investigations in Sect. IV Al however, showed that 
in general the situation is less obvious. For N ^ M 
Figs. El show a more complicated evolution pattern. 
Indeed, as we show next, in this case a new, shorter time 
scale enters the evolution. 

First, let us consider the simplest example, an obvious 
limiting case N = 1. As can be easily seen from Eq. I|16JI . 
the evolution contains just a single oscillation frequency, 
as J min = M/2 - 1/2 and J max = M/2 + 1/2. The 
corresponding value of r](M/2 - 1/2) is 2/{M + l) 2 . We 
thus have 



Pi(t) = 1 



+ 



(1 + Af) 2 (1 + M) 



■cos[gt(M + 1)]. (50) 



We see that the oscillation timescale grows linearly with 
the number of spins, not as a square root as in Eq. J3HJ|. 
We also notice that the amplitude of the oscillation is 
suppressed by the second power of the number of spins. 

Next, let us consider a more nontrivial case y/N + M < 
\N - M\ < (N + M). This case is characterized by the 
decay of the wavetrain to the value Pi given in Eq. (|4fcj|) . 
It turns out that for \JN + M < \N—M\ <C (N + M) the 
decay behavior is qualitatively different from Eq. I|48|l . 
The wavetrain in this case is approximately described by 



Pi(t) 



1 - - 




(2gtA 



2 arctan 



g 2 t 2 YP- 



(51) 



For details see Appendix iDl where an asymptotic series 
for Pi(i) is derived. Eq. I|51|l contains the leading term 
of this series. For sufficiently large \M — N\ (compared 
to vN + M) it gives an excellent approximation of the 
true result, as can be seen in Fig.[S] 
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FIG. 6: Time evolution of Pi(t) in the case N = 2300, M = 
2700, mjv = N/2, m M = -M/2. The Figure illustrates that 
Eq. 1511 (curve) provides a very good approximation to the 
exact solution obtained by summing the series in Eqs. 1161 . 
1A3H . 1A4H (dots). The time variable is defined as r = gt(M + 
N), the same way as in Figs. |H] El 

By inspecting Eq. (|51|) , we clearly see that for nonzero 
\M — N\ a new frequency, 2gA = g(N — M) enters 
the problem. We stress that this frequency depends on 
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N — M linearly. The amplitude of the oscillations at the 
maximum is suppressed by (N + M)/\N — M\ 2 , so as a 
function of \M— N\ the oscillations increase in frequency 
and decrease in amplitude. Finally, the oscillations de- 
cay with time as only a power law and the decay time, 
t eq ~ \N—M\/(g(N+M)) increases with \M-N\. Com- 
pare this with the situation for N = M, in which case 
the decay time scale is t eq ~ l/(gy/~N) and the decay has 
exponential dependence. 

At last, what conclusions can be drawn about the ap- 
plicability of the one-particle coherent Hamiltonian in 
this case? Our conclusions is that one can define in what 
sense it is applicable, but the argument is a bit subtle. 
Unlike in the case N — M considered in F&L II, in the 
more general setup considered here we find a new fre- 
quency, linearly dependent on the difference \N — M\. 
Can this be reconciled with Eq. (|15fl . which states that 
in the neutrino system which begins as a collection of 
flavor eigenstates there should not be any flavor coherent 
conversion? We recall that while the frequency of our 
solution increases with \N — M\, the amplitude of the 
corresponding oscillations becomes smaller and smaller 
(see Eq. H511 ). The oscillations become high-frequency as 
the system enters a state of freeze-out. The one-particle 
coherent limit corresponding to neglecting the residual 
small oscillations. 

A different way to state this is that the oscillation am- 
plitude is suppressed by (M + N) / (N — M) 2 and becomes 
large only when \N — M\ approaches N + M. In the 
latter case, however, the oscillation frequency, obviously, 
becomes of order \fN + M, and hence indistinguishable 
from the incoherent time scale. 



C. The Period of the Probability 

In this section we study the periodicity of the proba- 
bility P\(t) observed in Fig. EI Recall, the probability is 



JJV+J2-1 

P(t) = const + ^2 g{ J) cos [gt (2 J + 2)}. (52) 

J=\3n-3m\ 

The periodicity is a consequence of the probability being 
a sum of cosines with frequencies that are multiples of 
the lowest, fundamental frequency. 



1. The Case for N=l or M=l 

We first consider the case when N or M is one. Con- 
sider N = 1. For this case there is only one cosine term 
with J = ^ - |. Therefore Eq. reduces to, 

jN+j2 fM-l\ 
P(t) = E °( J ) + »? ( ~ S~ ) cos & ( M + !)] 



Hence, 



J=\3N-3m\ 



V 2 



T N= 



2tt 



g(M + l) 



(54) 



If we take M — 1 instead, the result is the same except 
N and M are swapped. 



2. The Case for N > 1 and M > 1 
Each cosine in Ea. l|52(l satisfies, 



cos [gt (2 J + 2)} = cos [g(t + Tj){2J + 2)} 



(55) 



Here Tj is the period of the cosine corresponding to 
angular momenta J. Now for all J, 



Tj = 



2tt 



g(2J + 2) 



(56) 




g(N+M) ' 



(53) 



To find the period we need the least common multi- 
ple of the Xj's. We find the period (together with the 
previous subsection) to be, 



if N + M = even (or J=integer) 
if N + M = odd (or J= half integer) 
if N=l or M=l 

(57) 

Note that the period is the same for the case of a flavor 
eigenstate background and a flavor superposition back- 
ground because the period only depends on the argument 
of the cosine. The argument of the cosine is the same for 
both cases. The discontinuity between the period when 
N = 1 or M = 1 and the other cases arises because in 
the first two cases there is an interference of many cosine 
waves (as there the cosines are summed over) and in the 
last case there is only one cosine wave. For M — N, our 
results reduce to that found in F&L II. 

Notice that, up to the factor of two controlled by 
whether N + M is even or odd, the period of the spin sys- 
tem depends only on the value of the spin-spin coupling 
g and not on the size of the spin system. For a neutrino 
system, g — v y' F and the period seems to depend on 
the volume occupied by the neutrino system. In fact, in 
a real neutrino gas, non- forward scattering effects would 
destroy any periodicity. We recall that for sufficiently 
large numbers of neutrinos the time scale of incoherent 
forward scattering, Eq. I|14(l , is much smaller than the pe- 
riod found here and, moreover, that forward incoherent 
scattering is only a small fraction of the total incoherent 
scattering (see F&L II). Thus, the periodicity is an ex- 
ample of an effect in the spin system that is not realized 
in the actual neutrino system. In the context of the spin 
system, on the other hand, the effect is perfectly physical 
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and indeed is reminiscent of the well-known "spin echo" 
effect experimentally observed in actual spin systems. 

As a final comment note that when we approximate 
the sums over the cosines by the integrals, in Eqs. 1)41)0. 
|@BJ, and (|5*TJl . the periodicity of the solution is destroyed. 
This has to do with the disappearance of the fundamental 
(lowest) frequency in the system. The same observation 
was made in F&L II for N = M. 



3. A note about minima 

Let us call the case when Pi = 1/2 + 5ZjC(J) — 
J2jV(J) (which is the lowest possible value that the 
probability can be) a perfect minimum. For the perfect 
minimum to occur we must have, cos [gt (2 J + 2)] = — 1 
simultaneously for all J. The period of each cosine de- 
creases as J increases, therefore we only need to find the 
time when the cosine with the largest period and the co- 
sine with the smallest period are simultaneously equal to 
-1. 

By a calculation analogous to that for the period one 
can show that if N + M = even there will never be a 
perfect minimum. If N + M is odd the times when the 
probability attains a perfect minimum is, 



T m in — — where fc = l,2,3... (58) 
9 

Note that this is half way between perfect maxima (the 
case where P r (t) = 1/2 + £j C(J) + T,jV(J))' Hence 
if M + N is odd (or J is half an integer) the probability 
has both a perfect minimum and a perfect maximum. 
This result shows yet another intriguing physical feature 
of the spin system. While a system with N+M = even is 
characterized by a set of recurring maxima, in the system 
with N+M = odd every other such maximum is replaced 
by a minimum. This behavior is illustrated in Fig. \7\ 
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FIG. 7: Plots of TV = 2560 and various numbers of spin down 
particles, M. The time is scaled so that r = gt(M + N). Note 
that for all the graphs N + M =odd 



Note that in this figure, as before, the abscissa is the 
scaled time t = gt(M + N). 



VI. THE FREEZE-OUT EFFECT AND A REAL 
NEUTRINO GAS: A CRITICAL ANALYSIS 



A. Overview of the problem 

In Sect.Elwe found that a general U N + M" (N ^ M) 
spin system that starts polarized in the z direction typi- 
cally equilibrate very little, even on the longer incoherent 
time scale, except in cases when /V and M are very close 
(\N — M\ < yN). Obviously, it is important to un- 
derstand if this finding corresponds to the behavior of a 
real neutrino gas. This question is not a trivial one. As 
mentioned in the introduction, going from a neutrino gas 
to the spin model involves several simplifying assump- 
tions: only forward scattering is kept, equal interaction 
strengths are assumed, the interactions are taken to be 
continuous in time, etc. All these assumptions could, in 
principle, introduce certain artifacts in the system and 
it is not a priori obvious that the freeze-out behavior 
found in the toy model is not just such an artifact. Wc 
will examine this question next. 



B. Interference of neutrino scattering amplitudes 

Consider a test neutrino, taken for definitcness to be 
\v e >, flying though a box with the gas of neutrinos. As- 
sume the neutrinos in the gas are in the flavor eigenstates, 



|5(0)) = \ v e v e ...v e v ll v IJi ...v l j) . (59) 

(N-l) M 

Let us limit scattering to only the forward process and 
ask: what is the probability that scattering changes flavor 
of the test neutrino? 

Let us review the arguments of F&L I. Omitting the 
flavor-preserving neutral current piece, which is the same 
for all neutrinos, we can write the result of the interaction 
to first order in the interaction amplitude a as 
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\v f ,)\v e v e ,..v f ,v ll v IJ ,...v ll ) — ► 

(N-l) M 

\F) = (1 + (N - l)ia)\v e )\v e v e ...v e v ll y^...v ll ) 

+ ia(\v ll )\v ( ,v e ...v e v e v ll ...v^) + \v^)\v e v e ...v e y ii v e ...v t i) + ...), 



(60) 



A I 



The interaction amplitude a is given by a = CSt = 
-V2G F {1 - cos@)5t/V 0, E3, where St is the dura- 
tion of the interaction between a pair of neutrinos. Since 
the final states are mutually orthogonal, the probability 
of flavor change goes as Ma 2 indicating the incoherent 
nature of the evolution. 

Eq. I|6l)[l assumes that the time step St is small such 
that not only CSt <C 1, but also NCSt <C 1 are satisfied. 
While the former is indeed true for any reasonable neu- 
trino system of interest 01 > the latter need not hold, as 
Na or Ma for large N and M are not necessarily small. 
To make contact with the spin system, which stays con- 
tinuously coupled, we must consider longer interaction 
times. Correspondingly, we need to modify Eq. (|60ll . 
keeping higher powers of a. 

One can accomplish this by regarding (|6U[) as the first 
order term in the perturbative expansion of the wavefunc- 
tion at the time St (NCSt < 1, MCSt < 1), and using 



this to make further steps in time to the time T = nSt 
[CT < 1, but not necessarily NCT < 1, MCT < 1). 
Let us define a convenient notation of the initial state 

1 1) = |t/ e )l lyy^ i/^),...^) (61) 

(N-l) M 

and the "exchange" state 

\E) = \E 1 ) + \E 2 ) + ..., (62) 
, ' 

\Ei) = \v li )\v e v e ...v e ,v e v lJL ...v li ), (63) 
\E 2 ) = \v li )\v e v e ...v & v lJL v e ...v li ), (64) 

Then at the time T the state is 



\T = n5t) = (l + i(JV- l)CSt) n \I) 



iCSt 



J2 (! + i(N - l)CTt)"- 1_r (1 + i(M - l)C8t) r 



0{[iCStf 



\E) 



(65) 



as can be readily verified by induction. 



The cyclotomic polynomial is summed to give 



\T = nSt) = (1 + i(N — l)CSt) n \I) 

l (1 + i(N - l)C5t) n - (1 + i(M - l)C5t) n 



iCSt- 



0{[iCSt} 2 ), 



iC(N - M)5t 



\E) 



(66) 



r 



and, taking the limit n 



\T = nSt) 



i{N-l)CT 



\I) 



,i(N-l)CT _ i(M~l)CT\ 



(N - M) 
0{[iC8t] 2 ). 



-\E) 



(67) 



Note that the explicit dependence on St has now disap- 
peared. This is to be expected since the formal solution 
of the Schrodinger equation for the system depends only 
on the time difference between the initial and final times. 
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The reason for the freeze-out behavior is now clear. For 
large values on N — M, the probability of a transition to 
the exchange state is suppressed by the factor (N—M)~ 2 . 

Notice first of all that this effect disappears when 
N = M. This is why in this case one can observe full 
equilibration and test the incoherent nature of the evo- 
lution. However, one should note that in the present 
approximation the N — » M limit gives iCT as the am- 
plitude for the state \E). As in this case the phase of 
the initial and exchange states are the same, we regard 
one factor of T in the probability as a relic of the energy 
conservation delta function, and obtain as the probabil- 
ity per unit time for a transition from the initial state 



§ = M|C|2 



(68) 



(remember that (E\E) = M), showing incoherent equili- 
bration. 



C. Does the freeze-out happen in a real neutrino 
system? 

Finally, we can address the key question: does the 
freeze-out observed in the model of interacting spins ac- 
tually happen in a real neutrino system? The answer is 
negative. 

Let us return to the arguments in the previous sub- 
section. Upon closer inspection of Eq. H6U|) . it becomes 
clear that for large N, M and times T <C C _1 , the prob- 
lem is in fact equivalent to an oscillation set-up with 
a Hamiltonian that, in the basis \I), \Ei), \E 2 ), has 
(iCN,iCM,iCM,iCM, ...) on the diagonal 2 and factors 
of iC on the off-diagonal, in the first column (and the cor- 
responding terms in the first row). For the case N = M 
considered in the transition probability between |J) 
and any of the states \Ei) is given by P ~ \itC\ 2 . The 
generalization of this to N ^ M is 



ixC x / dt exp[i(7V - M)t x C] 



(69) 



It is straightforward to check that the transition prob- 
ability that follows from Eq. (|69|l 



exp[iC(N - M)T] - 1 
(N - M) 
2 -2cos((JV - M)CT) 
(N - M) 2 



(70) 



agrees with that obtained from Eq. I|67l) . indicating that 
the corresponding amplitudes are the same, up to an 
overall phase convention. 



Let us take the time T in Eq. (|69|l to be the timescale 
of incoherent equilibration, given by Eq. I|14|) : 



(G F y/(N + M)/2/VY 



(71) 



If N = M, the probability to transition into any one of 
the exchange states \Ei) in this time is P ~ (Ct*" c ) 2 . Re- 
membering that C = — v2Gf /V(l— cos 9) and dropping 
the angular factor (1 — cos 8) (which is also dropped in 
the derivation of Eq. l|7T|)l. we get P ~ ((N + M)/2) -1 . 
The probability of transition into any of the exchange 
states is P(\I) -> \E any )) ~ (TV - 1 + M)P ~ 1/2, 
thus confirming that the equilibration does happen on 
the time scale tVl c . 

If, however, we have (N — M)CT > 1, the oscillat- 
ing exponent in Eq. (|69|l suppresses the transition and 
the system freezes out. Thus, we immediately obtain a 
physical estimate for the freeze-out: 



\N - M\y/2/(N + M) > 1, 



(72) 



We assume large N and M and neglect terms of order one. 



which is nothing but the condition we found in the spin 
system, Eq. (I45f) . up to a trivial factor. This provides a 
powerful check of the validity of Eq. (|69l) and the physical 
picture that led to it. 

We can now give interpret Eq. (|69|) the following inter- 
pretation. For N, M 3> 1, the states \I), \Ei) are approx- 
imately energy eigenstates. Strictly speaking, of course, 
the true energy eigenstates are mixtures of \I), \Ei), but 
for large N, M one of the eigenstates is predominately 
1 1), with a small admixture of the other states, and vice 
versa. Correspondingly, we can loosely speak of the en- 
ergies of the | J) and \Ei) states. For N ^ M the state \I) 
and any of the states \Ej) have different energies. The 
incoherent spin flips in the spin system are thus forbid- 
den by energy conservation (enforced by the oscillating 
exponent in Eq. H69|) ~). 

In a real neutrino system, of course, there is an addi- 
tional degree of freedom: the neutrino momentum. The 
energy liberated in the process of flavor exchange is con- 
verted into a slight shift of the neutrino kinetic energy 
and the energy conservation is assured. In the simplest 
manifestation of this point, a neutrino existing a dense 
neutrino region will gain a small amount of kinetic en- 
ergy. (The neutrino-neutrino interactions, just like any 
interactions between like particles mediated by a vector 
boson, are repulsive.) 

Thus, the freeze-out observed in the spin model is an 
artifact of the model, more specifically, a consequence of 
the dropped degrees of freedom, neutrino momenta. 

This discussion applies to the system of neutrino plane 
waves in a box, which was the system for which the spin 
model was created. It is instructive to instead consider 
a system of interacting neutrino wavepackets. In this 
system the interaction energy is localized to the regions 
occupied by the overlapping wavepackets, with the to- 
tal interaction energy (averaged over volume) being the 
same as in the case of plane waves. The statement of the 
changing kinetic energy, as the neutrino moves from a 
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more dense to a less dense region, has its analogue in the 
case of wavepackets as well: The kinetic energy of a given 
packet changes as it interacts with another packet. Upon 
spatial averaging this change, one gets the "plane waves" 
result. There are fewer interactions in less dense regions, 
hence the averaged kinetic energy there is higher. 

The advantage of thinking about wave packets is that 
in this case it is clear that the rate of incoherent scatter- 
ing cannot be suppressed: for sufficiently small wavepack- 
ets they interact with each other two-at-a-time, giving 
the usual scattering rate. The answer in the physical 
system cannot depend on the size of the wavepackets, 
since the interaction Hamiltonian is energy-independent. 



VII. CONCLUSIONS 

No problem in physics can be solved exactly. Often, 
however, a given physical problem could be mapped onto 
an idealized system, for which a complete solution exists. 
One then needs to critically examine which features of the 
solution carry through to the original physical system and 
which arise as artifacts of the simplifying assumptions 
made along the way. 

We have performed just such an analysis. Our goal 
was to investigate coherent effects in a gas of neutrinos. 
We have simplified the system by limiting the scattering 
to only forward direction and omitting the momentum 
degrees of freedom. With only the flavor degrees of free- 
dom remaining, the neutrino system could be mapped 
onto a system of coupled spins. We further simplified 
the problem by setting the coupling strengths to be equal, 
and considering a certain class of initial states. We have 
shown that the resulting spin model could be solved ex- 
actly. Our solution generalizes the analysis of F&L II 
0, where a particular initial state, the one with equal 
numbers of up and down spins, was considered. 

The solution proved to be very instructive. The spin 
system exhibited both coherent oscillations and incoher- 
ent decay, reproducing effects expected for a real neutrino 
gas. Several examples were considered in Sects. IIVI and 
fVl As these examples demonstrate, coherent effects in 
the neutrino system correspond to the behavior of the 
spin system in the classical limit; likewise, incoherent ef- 
fects in the former have their analogue in quantum effects 
in the latter. 

We have presented two different methods of obtaining 
the solution. The first construction separates the system 
in three parts: the test spin, the other spins that started 
in the same orientation and the remaining spins. The sec- 
ond method exploits the high degree of symmetry of our 
chosen initial states. Both methods have complementary 
advantages. The second of the methods provided crucial 
physical intuition for understanding the classical limit of 
the system and the interplay of quantum and classical ef- 
fects, as shown in Sects. llVl and IVAl At the same time, 
the first method is more general and can be used to gen- 
erate solutions for any initial state of the test spin (by 



simply changing the state of the test spin by adjusting 
the parameter A in Sect. 1111 C 31 ). 

Several interesting results were found. One such result 
involves the subtlety of identifying coherent effects in the 
neutrino system by taking a classical limit of the spin sys- 
tem. The special case considered in pj seemed to suggest 
that in the large N quantum effects separate from classi- 
cal ones necessarily by having longer time scales (scaling 
as a square root in the number of particles). What we 
found in Sect. IVBl is that some quantum effects instead 
decouple by having a vanishing amplitude, while their 
frequency scales as if they were classical effects. 

Another interesting finding of our analysis is that in 
certain spin systems the evolution does not lead to equili- 
bration, even on the incoherent time scales. For example, 
the system considered in Sect. [V] is "frozen-out" , unless 
the initial numbers of spins up and down are very close. 
By studying this example and comparing it to the case 
considered in Sect. IIVI we conclude the freeze-out hap- 
pens in the spin system for any mode where quantum 
effects are suppressed a large classical (coherent) effect. 
We have seen, e.g., in Sect. IV Al that equilibration is pre- 
cluded by the conservation of the classical angular mo- 
mentum. Only when the classical angular momentum is 
reduced to the level comparable to the quantum effects 
is the equilibration allowed to proceed. The special case 
considered in falls into the latter category. 

We have explained in Sect. I VII that the physical origin 
of the freeze-out effect comes from dropping the momen- 
tum degrees of freedom of the neutrino system. While 
the freeze-out certainly happens in the spin system, it 
does not happen in a real neutrino gas. This gives an 
important example of the limitations of the spin model: 
while it captures all coherent effects in the neutrino sys- 
tem, it does not always capture incoherent effects, only 
those that are not classically suppressed. 

In addition to giving the full solution in Eqs. l(TT)|) . l(T7|) . 

and 118|) , in several case we have derived approximate ex- 
pressions that elucidate the physics of the system. For 
example, Eq. I|40() cleanly demonstrates the presence of 
both oscillation and decay in the spin system. Eq. I|4t)|) 
helps to understand the interplay of classical and quan- 
tum effects in the processes of equilibration and freeze- 
out. Finally, Eq. I|51(l showed the subtlety of decoupling 
the quantum effect in the classical limit. 

We believe that our extensive investigation here will 
contribute to better understanding of the flavor evolution 
in dense neutrino systems. We also hope that our results 
could find applications beyond the neutrino field. 



Acknowledgments 

A. F. was supported by the Department of Energy, 
under contract W-7405-ENG-36. I.O. and B.McK. were 
supported in part by the Australian Research Council. 



18 



APPENDIX A: EXPLICIT EXPRESSIONS FOR 

C(J) AND r)(J) 

Eqs. JOU) C3>! CHJ give the general form of the so- 
lution for the probability P\{t). Here we give explicit 
expressions for C(J) and 17 (J) for the specific cases ana- 
lyzed in the text. 

In Sec. IIVI we consider the initial state containing N 
spin up particles (m^r = N/2) and M spins in the or- 
thogonal direction (mjvf = 0). In this case, Eqs. iJTTJl 
P&[l. upon evaluating the Clebsch-Gordan, 3j- and 6j- 
coefficients, become 



C(J) 



v(J) 



(1 + 2J)[A(1 + E) + J(1 + J)] 
4J(1 + J)r[l + J — AT/2]r[l + J + A 

r[i + j + A]r[i + j + iv/2]r[7v + 1] 



rp- j + i]rp + j + 2] 

r[2 + J + A]r[2 + J + N/2]T[N] 

(1 + j)r[i + j - 7V/2]r[i + J + A] 
1 



r[s- j]r[s + j + 2]' 



,(Ai) 



(A2) 



Here A 



_ JV-M 



and E = ^±M. 



Similarly, in Sec. [V] we consider the initial state con- 
taining N spin up particles (tojv = N/2) and M spin 
down particles (rriM = —M/2). In this case, Eqs. (fTTjl . 
(fT%)l evaluate to 



C(J) 



(1 + 2J)MA[A(1 + E) + J(l + ./)] 
2 J(l + J) 

r(M)r(jv) 
r[E- j + i]r[E + 7+2]' 
m[(i + j) 2 - A 2 )r(M)r(7V) 
(1 + J)r[E- j]f[E + j + 2] • 



(A3) 



(A4) 



Eq. (|A3|> is ambiguous for J = 0. In this case, which 
occurs only when N = M, we have 



(7(0) = 0. 



(A5) 



Formally, this can be seen from Eq. i|17|) . in which both 
the 3j- and the 6j-symbols contain objects that do not 
form closed triangles and hence vanish. 

Notice that our expressions agree with the solution 
given in F&L II [| in the case N = M. Indeed, set- 
ting M = N, x = 0, and y = N we have C(J) = and 
V (J) = (1 + J)N[{N - 1)\} 2 /[{N -J- 1)\{N + J+ 1)!]. 
Noticed that Eq. (4.6) in the journal version of [ij con- 
tains a typo: the factor (1 + 2 J) in the numerator should 
be (2 + 2J) (the first line of that equation is correct). 



APPENDIX B: DERIVATION OF EQ. dlHt 

The derivation of Eq. 14U|) proceeds as follows. Let us 
start with the definition of i](J) in Eq. I|18|) . Approx- 
imating the gamma functions by the Sterling formula, 



T(z) = (z-l)!~ 72^(2-1/2) in(*-i)-( 2 -i) 5 an d combin- 
ing the resulting exponents, we get 

V(J) = 7T-4wS= ex P[^]i 



2 + 2J + M-N 



(1 + J)V2^ 
A = i U3 + 2J + M- JV)ln 



+ (3 + 2J + A^)ln 

- (2J-iV + l)ln 

- (2J-M + iV + l)ln 
+ (2JV — 1) ln[AT — 1] 

+ (1 + 2J- M- N)ln 

- (3 + 2J + M + N)ln 

- 2). 



2 + 2J + N 



2 

2J- N 



2,1 - AI + N 



-2-2J + M + N 



2 + 2 J + M + N 



(Bl) 



We need to find the maximum of A and expand A in a 
series in J to the second order around this maximum. A 
derivative dA/dJ has two types of terms: logarithms and 
terms that go like l/#, where "#" denotes an expression 
linear in J, TV, and M. Clearly, for sufficiently large val- 
ues of N and M the logarithms dominate. The maximum 
is then found by combining the logarithms and setting the 
argument of the combined logarithm to 1. This yields 



(2 + 2J + M - N){2 + 2J + N) 
(2J- N)(2J - M + N) 
-2 - 2,1 + M + N 



(2 + 2 J + M + N) 



(B2) 



For large N and M, we can drop small numbers compared 
to J, M, N. The above equation then gives the answer 



VM 2 + TV 2 



(B3) 



which is exactly the answer for the classical problem. 

Upon evaluating the second derivative, d 2 A/dJ 2 , we 
again find two types of terms, some that go like l/# and 
some that go like l/# 2 . Keeping only the terms of the 
first kind, substituting Eq. I|B3|) . and simplifying, we get 



(B4) 



Finally, we need to find the value of A in Eq. (|BL|> 
for J = Jmax- This involves a rather lengthy calculation 
involving a lot of cancellations. All terms of the type 
Jln[#], AHn[#], and Mln[#] cancel out, when Eq. |Bl|l 
is used. The terms of the type ln[#] give ln(M 2 /(47V)) 
upon simplification. Finally, the order one terms also 
cancel out. Thus, we get 



= In 



AN + \N' M' J 



(B5) 
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>r VM 2 +n ,j ^ p or j ar g e ^ Gaussian is sufficiently 

2 narrow and we can replace the prefactor 1/(1 + J) by 

A ~ 1 In ( — J - ( — + ^ ) ( J - \/M 2 + A 2 /2) 2 2/VA/ 2 + A 2 . The sum over J in Eq. © can be re- 

2 \4A/ \A M 2 y placed by an integral, in which we can extend the limits 

(B6) 

of integration to ±oo. We have 

This means that r](J) is a Gaussian centered at J ma x = 



/oo 
dJri(J) cos(2. 9 Jt) = 
-OO 



7lf2\ 1 / 2 /*00 (2 , 2N\( r VM 2 + N 2 \' 

' ' 1 dJe U + ^V 2 ) 



2ttVM 2 + N 2 \4A, 



cos(2g Jt) 



y/Hy/M 2 + TV 2 \4A 
1 M 2 



^2(M 2 + A 2 ) 



i(gVM 2 + NH) = 



2 (M 2 + TV 2 ) 



/ 2 i 27V \ — 1 2.2 / 

e -U+i^J 9 * cos(g\/M 2 + N 2 t). 



r 



(B7) 



The constant term j°j min C*( J) can be immedi- 

ately found as the difference between 1 and the oscillating 
term at t = 0: 



1 - 



M 2 



1 1 TV 2 

2 + 2 (Af 2 + TV 2 ) ' 



(B8) 



2(M 2 + A 2 ) 
This concludes the derivation of Eq. (|4(J|) . 

APPENDIX C: DERIVATION OF EQ. l|46|> 

To find the constant value Pi to which P\(t) relaxes 
we can either compute the sum 1/2 + C(J), or 



use the trick at the end of Appendix [B] and compute 
1 — Sj^T^ 1 v(J)- Let us do the latter. 

The first step is to approximate the gamma functions 
in Eq. l|A4fl by the Sterling formula and expand the ex- 
ponent. After a fairly straightforward calculation, one 
finds a Gaussian in J + 1 : 



T(M + l)T(N) 



T[(M + N)/2 - J]T[(M + N)/2 + J + 2] 



exp 



B 



2{M + N-1) 
(M + TV) 2 



where 



B 



1 



(CI) 



M + - logM + TV - - ln(7V - 1 



(1 + M + A)ln 



M + N 



(C2) 



For definitiveness, let us take M > N and let us assume 
that [M - TV) < TV. We can then expand Eq. (fUgJ) in 



series in A: 
B ~ 



.Y-i)lnLY-1i - { V + ^ ) log A' 



2 

2 TV - I 

8A 2 



(A ~ Mf 



-In A + —{N - Mf + <3(1/A 2 ). (C3) 



We have 



[(1 + J) 2 - A 2 ] 
(1 + J) (A + 53) 



exp 



A 2 



(C4) 



where, as before, A 



_ N — M 



andS= 



Let us approximate the sum over J by an integral 

(M+iV)/2-l 

A = i - X] *?( J ) 

J=(M-A')/2 



1 



*-\ T (J 2 -A 2 ) 

a J — — r exp 

J(A + E) 



A 2 



A 



(C5) 



Let us shift the integration variable J = J— A and extend 
the upper limit of integration to infinity: 



1-Pl 



- (J + 2A)J 

a J — exp 

(J + A)(A + E) 



(C6) 
(J + A) 2 " 



(A + E) 



Introducing a new variable £ = J/\/E, we get 



1-Pi 



exp 



-A 3 



(A + E)E 



(C7) 
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where 



F(a) 



1 (C + 2oK i 



-(£ 2 + 2£a)] (C8) 



The function F(a) drops to zero quickly as a is in- 
creased beyond one. Since the argument of F in Eq. (|C7|) 
is A/ yS, we conclude that indeed Pi — > 1 for A > -\/E ■ 

For A < VE the argument of the exponent in Eq. (|C7|) 
is of the order A 3 /E 2 < E -1 / 2 and hence in the limit of 
a large number of spins the exponent can be set to one. 

The integral defining the function F(a) can be further 
transformed as follows: 



F(a) 



d( 



exp [-(£ + a) 2 ] e a 



« + «) 



r 2i 

y exp[a J 



dy 

y 



Finally, we arrive at 



Pi 



1 - 



-F 



A 

E 



where 



W 2 2 



dy- 



exp [-(£ + a) 2 ] 
-exp[-(£ + a) 2 ] 
(C9) 

(CIO) 

(Cll) 



The function F(a) equals 1/2 at a = (F(a) a =° 1/2+ 
l/2(21na + 7 )a 2 + l/2(21na + 7-l)a 4 + ..., 7 ~ 0.577216 
being the Euler constant). As a — > 00 -F(a) approaches 



(F(a) 



V2 + ...)■ 



APPENDIX D: DERIVATION OF EQ. (JSTJ) 

In Sect. El we consider the evolution of the spin sys- 
tem in which all spins are initially in flavor eigenstates. 
In this Appendix we derive an approximate expression 
describing this evolution in the regime when \J M + N < 
\M-N\<s:(M + N). 

The time dependent part of the evolution is given by 
J2j V(J) cos(2gi(J+l)). The coefficient r}{ J) depends on 
A = (N — M)/2 and can be usefully approximated for 
the case when A is small by Eq. (|C4ll of Appendix El 

As before, we approximate the sum by an integral 



(m+jv) /2-1 

2 r7(J)cos(2 5 t(J+l)) 

J=(M-N)/2 

(J_j-A_j 
dJ — 7-. r exp 

J(A + S) 



(Dl) 



cos(2gtJ), 



shift the variable, J to J = J — A, and extend the range 
of integration to infinity. We find 



(M+AO/2-1 

53 i](J) cos(2gt(J + 1)) ~ exp 

J={M-N)/2 



-A 3 



(A + E)S 



1 



(A + E) 



(J + 2A)J 

a J — = exp 

(J + A) 



' E 



2JA 



cos(2 5 i(J + A)). 



(D2) 



F&L II considered a limit A = 0. In that case the integral 
equals 

l-^VNgtcxp(~Ng 2 t 2 )eYZ(VNgt), (D3) 

from which Eq. I|48(l immediately follows. 

We notice that the values of J that contribute in the 
case A = are of the order VE, as the integral in 
Eq. ljD2|) is cut off by the term exp(— J 2 /E). If we take 
A > \/E, however, we find that the integral is first cut off 
by the exponential term exp(— 2JA/E), at J ~ E/2A. 
If A is sufficiently large, we may evaluate the integral in 
Eq. I|D2|) with the first term in the exponent dropped. 



The integral then becomes 



~ (J + 2A)J 
a J — — — — exp 



(J + A) 



2JA 



cos(2 ff i(J + A)).(D4) 



We next take advantage of the fact that integrals of the 
form J Q dx x n exp[— x] cos(ie + a) can be done and the 
answer has a simple and useful analytic form 



/•OO 

/ dxx n exp [—Ax] cos[-B(x + a)] = 
Jo 



{A 2 + B 2 ) 



■ cos 



Ba + (1 + n) arctan ( — 



.(D5) 
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Expanding = 2— ^- + ^ — ■^■ + ..., we can reduce (|D5fl . Evaluating this series, we obtain 

the integral in Eq. (|D4|I to a series of integrals of the form 



J 



/(*) 



~(J + 2A)J 
a J — — — — exp 



(J + A) 



2JA 



cos 



2giA + 2 arctan 



y- ( _ ir +i - A i-„/ 2<2 



n=2 



A 2 
E 2 



cos(2gt(j + A)) 

(f) 

(n+l)/2 

cos 2gtA 



fgt^X 

[l + n) arctan I I 



(D6) 



The form of Eq. (|D6|> suggests that the solution is os- 
cillating in time. It is instructive to understand how this 
oscillatory behavior disappears in the limit A 2 /E — > 0. 
First of all, notice that the summation in Eq. I|D6|) goes 
to some n m ax, not to infinity. Indeed, the presence of 
the factorial n\ indicates that we are dealing with an 
asymptotic series. For any given choice of A and £ there 
will be an optimal number of terms in the series, n max , 
that approximates the original integral best. The terms 
beyond n max grow in absolute value and the series di- 
verges. We easily estimate n max from the condition that 
the ratio of the two consecutive terms at t — be ~ 1, 
(n/2)A- 1 A- 1 E - 1, or 

?W ~ 2A 2 /E. (D7) 

Clearly, for A 2 < £ the expansion breaks down. 

Next, recall that in this limit the integral in Eq. I|D4|) 
is not valid anyway (the first term in the exponent in 
Eq. 1D2|) cannot be dropped) and one needs to consider 
Eq. JD2J. The answer to the latter in the limit A 2 /£ -> 



is provided by Eq. (|D3|I . indeed not showing any oscil- 
lations. The oscillations thus appear as \M — N\ is in- 
creased beyond y/M + N, when the series in Eq. I)D6|) 
starts providing a better and better description of the 
true answer. 

The final answer for the probability Pi (t) in the regime 
\/M + N < \M - N\ < (M + N) is given by 

Pi(t) = l-^2r,(J)+J2v(J)cos(2gt(J + l)) 
j j 

where I(t) is given by the series in Eq. ID6I Taking the 
leading term in the series and setting A = in the prefac- 
tor of Eq. (|D8() we arrive at Eq. I|51() . This approximation 
turns out to be quite accurate, as Fig- EI which shows the 
case N = 2300, M = 2700, illustrates. 
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